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Abstract 

A renormalization group study of a scalar theory coupled to gravity through 
a general functional dependence on the Ricci scalar in the action is discussed. 
A set of non-perturbative flow equations governing the evolution of the new 
interaction terms generated in both local potential and wavefunction renoma- 
lization is derived. It is shown for a specific model that these new terms play 
an important role in determining the scaling behavior of the system above the 
mass of the inflaton field. 
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I. INTRODUCTION 



It is well know that modifications of Einstein's theory of gravity are required in order to 
include quantum effects. Although a consistent fundamental theory of the spacetime near 
the Planck scale is still not available, one can consider General Relativity as a quantum 
effective field theory fjl^]. From this point of view one gives up the requirement of (per- 
turbative) renormalizability to consider a more general action that is a functional of any 
geometrical invariant 5i which can be constructed from the principle of general covariance, 
3? = {R, R a pR al3 , C a p^C a ^ h , ...}, and of the matter field and their derivatives 

S[4>,g] = J <fx^=g~C{% 0, V 2 3?, V 2 0, V 2j K, V 2j 0,...,) (1) 

This new theory is analogous to the Euler-Heisenberg lagrangian, low-energy effective theory 
of the more fundamental microscopic QED. Due to the smallness of the ratio between the 
Newton constant and the Fermi constant, this modification of Einstein's theory is experi- 
mentally indistinguishable from the standard theory at ordinary energy scales ||]. However 
at higher energy scales [|J it has significant physical implications, mainly for the inflation- 
ary scenarios ||, nowadays tested by measuring the temperature fluctuations of the cosmic 
microwave background radiation (CMBR) over large angular scales. 

Higher-derivative gravity theories also arise through the coupling between a quantized 
field and the classical background geometry |J. In this case one adds new terms to the 
standard Einstein-Hilbert lagrangian. Those terms are general functions of 3? and are needed 
to cancel the ultraviolet (UV) divergences at any given order in perturbation theory. This 
approach is useful to study the short-distance behavior of the model, but it loses its validity 
if we are interested in the low energy behavior. 

In fact, the standard perturbative approach based on the scaling property of the Green's 
functions under a rescaling of the metric handles only a finite number of operators, those 
which are important for the UV fixed point. In this way one cannot follow the evolution 
of the irrelevant, i.e. non-renormalizable operators. Although the standard classification 
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of the interactions is given accordingly to the behavior of the renormalization group (RG) 
transformation near a fixed point, we use a more general definition, calling a coupling irrel- 
evant, marginal or relevant if, in cut-off units, either it gets smaller, or it does not run, or 
it grows under a RG transformation. The irrelevant operators, even if they are not present 
in the bare lagrangian, mix their evolution with the renormalized trajectory of the relevant 
couplings. Although near the gaussian UV fixed point their running is suppressed, they can 
behave in a quite different manner in other scaling regions thus deviating the renormalized 
flow. 

It is therefore important to trace the evolution of all the couplings constants generated 
by the renormalization procedure in order to decide what can be neglected and what cannot. 
If for example a new scale other than the cut-off is present in the problem, the scaling laws 
change near that scale and the standard UV relevant interactions might not be sufficient 
to describe the physics at the crossover. This happens in the O(N) model where, in the 
non-symmetric phase, the presence of singularity in the beta functions for the longitudinal 
components indicates that new infrared (IR) relevant interactions, due to the Goldstone 
modes, appear after the condensate has formed ||. 

Actually we deal with a similar situation when we move up in energy, going from an effec- 
tive theory defined at some low energy scale, towards the high energy region. In fact integrat- 
ing out the heavy quarks contribution in the Standard Model generates non-renormalizable 
vertices that are suppressed at the weak scale || but that can show up when the heavy mass 
threshold is reached fllPf . 

Although there is no way to localize the energy stored in the gravitational field, a local 
definition of a "mass" scale can always be given, even if no matter field acts as a source in 
Einstein's equations, i.e. T^ v = 0. A possibility is to consider m G = yf^j where \l/ 2 is the 
coulombian component of the Weyl tensor. For instance, in spherical symmetry it coincides 
with the square of the Weyl tensor, and for a Schwarzschild spacetime ^2 = m/r 3 . Since 
tuq ~ \Raf3-y&\ one can think of as the local effective mass producing the curvature in 
that point. In this respect the gravitational field acts as an external field that influences the 
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quantum dynamics within scales Iq ~ rriQ 1 . 

On the other hand, during the slowroll phase of the inflationary era, the only relevant 
"mass" scale is provided by H = a/a, a(t) ~ a§t m is the scale factor, and the quantum 
evolution of the fluctuations is stopped after crossing the horizon. The subsequent classical 
behavior is therefore dependent on the hierarchy between these two scales, the mass of the 
inflaton field, mi, and the scalar curvature R. 

An important issue is related to the running of the conformal coupling constant £ which 
couples directly the inflaton with gravity through the ^R<p 2 term in the action. The value 
of this coupling is determined in principle by the spectral index of the CMBR temperature 
fluctuations n s = 1.17 ±0.31 ||11|| , depending on the model of inflation considered. Although 
for some models negative values of £ are preferred |12[ , in perturbation theory with the min- 



imal subtraction scheme one finds £ = 1/6 to be an infrared fixed point for the renormalized 
coupling |7|,|l3|. 



The renormalization group approach used in Statistical Mechanics is the best tool to 
study problems where many scales are coupled together WM. In this paper we shall use 



the differential form of the RG transformation formulated by Wegner and Houghton [|T5 
Starting from a bare action Sk at the cut-off k one first calculates Sk-Ak m 



e 



J?*-AkM = / V[i)}e- Sk[(i,+ ^ (2) 



by using the loop expansion, where ip and (f) respectively have non-zero Fourier components 
only in the momentum shells k — Ak < p < k and p < k — Ak. The differential RG 
transformation is obtained by taking the limit of an infinitesimal shell Ak/k — > bk/k. The 
higher loop contributions in Eq. (0) are suppressed as powers of Ak/k in the limit Ak — > 
for finite k and an exact, non-perturbative, one-loop RG equation is obtained 

k ^k~-~2 [[n ^^ ) + { ^y^^) ^ } {6) 

where the brackets indicates sum over the Fourier components within the shell. This func- 
tional equation rules the evolution of all the interaction terms that are generated in the 



renormalization procedure. Despite the solution is not know in the general case one can de- 
rive approximate, non-perturbative evolution equations in terms of the gradient expansion 
fl6|j , by writing 

S[<j>]= fd D xJ2U n (<P,d 2n ( p) (4) 

where U n is an homogeneous function of the field and of its derivatives. In this way Eq. ([3]) 
decouples in a set of infinite non-linear partial differential equation for the U n . The explicit 
construction of the RG transformation can be achieved by the blocking transformation or 
"coarse-graining" of the fields by means of a smearing function that introduces a sharp cut- 



off in the high momentum modes |17 |. A perturbative construction of the coarse-graining 
approach has been firstly discussed in a gravitational context in |L8| and proposed in the 
formalism of the "average action" in |T9|j20|| . In particular in ||19|| , by considering the con- 



tribution of the blocked potential Uo(4>) in Eq. (|]), it has been shown that £ = 1/6 is not 
necessarely an infrared fixed point for the massive theory coupled with gravity. 

An equivalent construction of the exact RG equations in the Wilsonian approach is 
presented in [|21|| . In particular in [^] the infrared behavior of the Einstein-Hilbert action 
is studied in the case of pure gravity. 

In this paper we explicitly construct the RG transformation for the blocked potential 
Uk and the wavefunction renormalization Zj, which are general functions of both <ft and the 
curvature scalar R. In particular we construct the RG traformation by means of an 0(D) 
symmetric smearing function (D is the dimension of the spacetime). This will be achieved 



by introducing locally a momentum space [23| and by working up to first non-trivial order 



in the expansion of the metric. The set of equations we derive generates the running of any 
interaction term of the form gijE^cfP in both Uk and Zk- 

We study the RG flow in the two different scaling regions of the ultraviolet and infrared 
domain. In particular the crossover in corrispondence of the mass gap is analyzed in a 
specific model by means of a numerical investigation of the equations. 

The organization of the paper is the following. In Sec. II the flow equations for Zk and Uk 



are derived. In Sec. Ill we analitically study the behavior in the UV and in the IR region. In 
sec. IV we apply the method to a model by numerically integrating the flow equations. The 
results of the numerical analysis are then displayed. Sec. V is devoted to the conclusions. 



II. RENORMALIZATION GROUP TRANSFORMATION 

In the following we suppose that the spacetime has a well defined Euclidean section Q. 
The general form of the Euclidean bare action, quadratic in the derivatives of the scalar 
field and of the Ricci scalar R is 

S A [<f>,g) = J d D x^{^Z^,R)d^d^ + V( ( p,R) + 

\w{<f>, R)d l j i Rd >1 R + R)d^R) (5) 

where A is the ultraviolet cut-off of the theory, V is the local potential, while Z, Y are general 
functions of and R. The action (||) is equivalent to the following one 

S A [<f>,g] = J d D x^g{~Z{ ( j ) ,R) ( j>V 2 cj> + V{ ( j ) ,R)- 

^W(<P,R)RV 2 R + y(<P,R)d^d»R} (6) 



provided 



where V 2 = V M V M is the Laplace-Beltrami operator, and the prime indicates derivation 
with respects to (p. 

In order to construct the renormalization group transformation in the Wilson-Kadanoff 
approach one first defines the average blocked field by coarse-graining the original field 



4> k (x) = J d x'y/gp k {x,x')4>(x') (8) 

where p k {x, x') is a smearing function which is constant in a given volume where the field is 
averaged, and it is zero outside. The blocked action is by definition 



e -s k [*, 9 ] = I D ^ l[5((f) k (x) - $(x))e- SA ^l. (9) 

"* x 

In this way one "averages out" the high energy degrees of fredoom that are relevant for the 
short-distance behavior of the theory. 

Let us consider the decomposition of a generic function f(x) of £ 2 (Q), in eigenfunctions 
of the Laplace-Beltrami operator ipi(x), f(x) = J2i fi^iix), where the ipi(x) represent an 
orthonormal base in this space (i is an integer if the manifold is compact). Once the 
spectrum of the operator V 2 is known, the sharp momentum smearing function defined in 
fjl| can then be used to explicitly construct the projector Pk(x, x') in Eq. (H). In particular, 
for a general spacetime one can introduce locally a momentum space and use the 0(D) 



sharp-cutoff smearing function as in |L7|. This is a good approximation as long as the cut- 
off of the theory is much above the scale of energy where the curvature becomes dynamically 
relevant and, at the same time, it gives a precise meaning to the UV cutoff A introduced 
above. The field is then split 

4>{X) = + £{X) = J2 #mll>m{x) + £ ^n{x) (10) 

M N 

where M and N respectively indicate the slow and the fast variables; more precisely M is 
the set of index values m labelling components with low momentum: p < k and N is the 
set of index values n and n' corresponding to high momentum k < p < A, k being a fixed 
momentum scale below the cutoff A. The functional integration in Eq. (§) yields 



S k [*] = -In / 2W>fc]n*(&(*0 -Hx))e' s ^ = 

J X 

-In [ Vi$]l[S($ m - $ m ) / V[£]exp{-S[#] - -W^K^ ...} = 

-In / mUWm - $ m )exp{-5[^] + WFnK^ - W (lniT) njn } = 

J M N Z N 

Sm -J2F n K-^F n/ # + ^Tr \nK[*) (11) 



N 

where 



K 



5 2 S 



(12) 



In coordinate representation we get 



e q -i 

F(x) = -— = — Z' df&dPQ - Z V 2 $ + V (13) 



and 



6<p(x) 



r2 c 

K(x,x') ' ' 
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8(t){x)8(t){x') 

{- x -z" dfidP® - z' {dpQdP + v 2 $) - z v 2 + v"}5{x, x') = 

V 2 $ - Z V 2 + V"}5(x, x') (14) 



where 



5(x, x') = g- 1/ \x)5{x - x')g~ l/ \x') = 5{x - x')/jgjx) (15) 



In deriving Eqs. (0,0), we have supposed that, for our purposes the fluctuations of R are 
negligible or, in other words, we have suppressed in the action the terms containing d^R; 
as a consequence W and Y, although in principle could be not vanishing, do not appear in 
Eqs. ( |Il3||l4"D and in the following calculations. 



The last equality in Eq. ( 13 ) follows after some manipulations, but the easiest way to 



verify it, is to compute the second functional derivative from the action (|(J), and use the 
relation for Z in Eq. (|7|). 

In order to explicitly carry out the blocking procedure, the functional form of the blocked 
action in Eq. (|TT|) is needed, and the usual assumption is to mantain the same form of the 
starting action; thus, according to the restriction to constant R, the blocked action reads 

S k [$] = J d D x^g^Z k (<S>, R) d^d^ + f4($, R) + 0(V 4 )^) (16) 

The function Z is generated by performing the path integral in Eq. ( |11~D for non constant 
field configurations 

$(x) = $ + 4>{x) (17) 
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By expanding the action in Eq. (|16|) up to second order in 0, one can single out the 
contribution to the blocked potential U and to the wavefunction renormalization Z, as it 
has been shown in [plj and in || . In a general spacetime this procedure might be problematic 



because the effective potential is in general ill defined when the spacetime dynamics becomes 
relevant. For this purpose it is convenient to consider energy scales much higher than the 
characteristic curvature scale k 2 <C R, and as anticipated above, to introduce locally a 
momentum space |2^] and to use an 0(D) symmetric smearing function in Eq. (H) to 
evaluate Eq. ([11]). In fact, although several techniques are available in the literature |24| to 



compute the functional determinant in Eq. (0), they are often combined with dimensional 
regularization and minimal subtraction prescription. Instead, the explicit use of a sharp 
cut-off regulator allows a clear separation between the UV ad the IR domain of the theory. 

Let us briefly outline the computation of the blocked action. It is convenient at this 
point to introduce a covariant notation: for any operator A the Lorentz invariant operator 
A is defined by A(x, x') = g 1 ^ i (x)A(x, x')g 1 ^ 4 (x'). In particular for the trace in Eq. fllTD, 
properly written in coordinate representation, we get 

Tr In K = J d D x^g~{x) hi K(x, x) = J d D x In K(x, x) = T^K (18) 

where the definition of the logarithm as power series of operators has been used. 
Let us examine the various pieces entering K. According to the above definition 
g 1 ^(x)V 2 x 8(x 1 x')g 1 ^(x') = V 2 x 5(x — x') where Eq. (j!5|) has been used and 

V 2 =g^(x)dl u + b^x)d, + a(x) 

&"(*) = g- 1/2 (d^) + 2g l l\d»g- 1 ' i ) + {dyfT) 

a(x) = g- x l\d ll Jg){&>g- 1 ' 4 ') + g^d^g' 1 ^. (19) 

Analogously 

g 1/ \x)Z\V 2 x ^)5{x, x')g 1/4 (x') = g- 1/ \x)Z' V 2 x (g 1/4 $)8(x - x') (20) 
Thus, making the replacements 



$(x)g 1/4 (x) -> ^(x)# 1/2 (x) -> V(x) (21) 

we finally get 

K(x, x') = {- l -Z' V 2 ($) - ZV 2 + V"}5(x - x') (22) 

where the prime now indicates the derivative with respect to the new scalar field of Eq. (pip . 

Note that we could have obtained Eq. (|22|), making the replacements of Eq. (|2~1|) in the 
original action and taking second functional derivatives with respect to the new scalar field. 
Since the following developments are naturally expressed in terms of the replaced variables 
in Eq. (pT|) , we shall keep on using the new field and potential, but for simplicity we shall 
not change the notation. 

We can now introduce normal coordinates {y a } around the point x' corresponding to 
y a = 0, expand the expressions in Eqs. (|T9|) and (|22|) up to some order and invert the 
propagator as shown in [2j| and [2l| where they fixed Z = 1 at the bare level. However for 



our purposes it is more convenient to work out separately the expansions for the metric and 
the field $ since in principle they involve different scales. We expand the metric tensor up 
to second order in its derivatives therefore writing 

g^( y ) = et + l -R\ v py a y p + 0(d 3 g) 

g = l- l -R aP y a y p + 0(d*g) (23) 

This truncation is consistent with the choice of neglecting d^R, that would appear starting 
from the third order in the metric expansion. 

In order to compute Z we first expand Eq. (|22|) around $(y) = $o + 4>{y) up to second 
order in obtaining K = K$ + K-^ + K^ 2 where (Zq and Vq are used for constant scalar 
field $ and □ = 5^8^ ) 

K {y,0) = {-Z D + - ^(2Z + $o4)}%) 
K^y, 0) = {~Z' (nj + 20Q) + - ^(3^o + *o3?)?}*(v) 
Kj 2 (y, 0) = {~ZZ(M + 2 n) + Ic> 2 - + «")0 2 }<%) (24) 
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In deriving these expressions we have taken into account that V , when applied to a scalar 
quantity, becomes V 2 = (□ + (1/3) ^//^V^ ~ (l/3)i^ a y a <9 M + ±R) = □ + R/6. We 



rewrite Eq. (18) using 



oo I -r \n+l 

Tt\nK = Tt\nK + j: { -^— Tr{K n (K^+K^) } (25) 



n=l ™ 



One therefore also expands the expression ( fCB]) up to second order in cf> and compare it with 
Eq. (PD, obtaining 



S k {§} = sl + s k + si = s A + s A + si + 

±Tr\nK + ±Tr K ~%+ \^ K^K^ - -^{K^K^K^ + O(0 3 ) (26) 

where S^S^S* 2 are, respectively, the zeroth, the first and the second order terms in the 
functional expansion of the action in The zeroth order gives the blocked potential 

Sl = Sl[%] + l -Ti\nK (27) 

wherefrom (in the following we use the notation J p = f k< \ p \ <A (f^p) 

Uk = Vo + ll p ln {Zop2 + V ° ~ l2 {2Zo + ®° Z ' o)) - (28) 



In deriving the last relation we have evaluated the trace in Eq. ( [27D in the Riemann frame. 
Since the integrand is calculated in $ an d K is diagonal in the local momentum represen- 
tation, a volume term can be factored out, yielding a covariant equality between constant 
quantities. The blocked potential U ($) is obtained by replacing the constant field $o with 
$(x). It should also be observed that Eq. (|28|) for k = gives the effective potential, ob- 
tained by taking the coincidence limit in the fluctuation determinant and by retaining only 
the first order terms in the Schwinger-De Witt proper time expansion [ ffi5[| . 

A similar procedure must be used for selecting the contribution of Z k in the other terms 



in the expansion (p6|). In this case, for computing the traces, according to |16| we retain at 
the same time operators in the momentum and in the coordinates representation, as long 
as the operators in one representation are on the right side of the operators in the other 
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representation. In the terms where the required ordering is not fulfilled, we employ commu- 
tation rules for p-dependent and x-dependent terms in order to move the p-dependence on 
the left of the x-dependece and disentangle the volume integration from the p integration. 
The main commutation rule, from which all the others can be deduced with no ambiguity, 
is 

[$(x), Pll ] = -td^(x) (29) 

Thus we get 

T7 Ko~% = J d D y J p ±(Z>p 2 + Vf - ^(3Z> + %Z^(y) 
TtK^K-^ = Jd D yJ ±(l Z '>{d^{y)d»4>{y) +p 2 ft(y)) + 



\<'9(y) - + «")0 2 (z/)) (30) 



where 



K p = Z oP 2 + V£ - (R/12)(2Z + $ ^o) (31) 

and, after having used the commutation relation (f29"D for computing Tr ln(Ko 1 K-^K Q 1 K~) 
we eventually obtain the following renormalization group equations 



U k = V+l { \n(K p ) 
2 Jp 



9 7 lf ( Z" 2Z'K' p Z{K' v f ZZ'p 2 K' p ZY {K' p )\ 

Note that there is no contribution from the term FK~ 1 F in Eq. (|TTD in the gradient 
expansion as it can proved by direct calculation. The reason is that in this term no volume 
integration occurs and one obtains the fourier transform of <p which is constrained to have 
p < k. Since the momentum integration is performed for p > k this term identically vanishes 

These equations have been obtained in the independent mode approximation where one 
computes the blocked action by an independent integration of the degrees of freedom between 
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the UV cut-off A and k. In the limit k — » one recovers the effective action. In particular 
if we neglect the field dependence in Z, therefore reducing it to a running coupling 77, the 
previous equations read 

u>= v+ u*y +v "- wm 



2J P (i)p 2 + V")|j,«,o 

Z ' = " + Uhg ^— )■ (33) 

where now if p = 77/j 2 + V^" — (R/6)r] and the cosmological constant that corresponds to the 
constant part of V has been modified in order to include the field independent part in the 
logarithm in Eq. fl3"3|). 

As an example let us consider the standard A0 4 bare scalar self-interacting theory in 
curved spacetime 

C = -rjd^d 1 *® + im 2 $ 2 + -£R<& 2 + ^A$ 4 (34) 

the RG equations (j3"3|) for Z can be integrated up to k = and one obtains in D = 4 the 
finite contribution to the wavefunction renormalization function 

1 A 2 $ 2 

Zk=0 = V + 1927T 2 m 2 + § $ 2 + g - p (35) 



which coincides with the result in 26 



We want to remark that, due to the replacements in Eq. the couplings and the field 
appearing in the bare lagrangian do not coincide with the ones appearing in the evolution 
equations, although, in case of small curvature the difference becomes practically negligible. 
However it is formally correct to expand in terms of running couplings the modified potential 
defined in Eq.(pTl) and not the bare lagrangian as in Eq. (3~1). 



III. NON-PERTURBATIVE FLOW EQUATIONS 



The previous equations have been obtained by eliminating the degrees of freedom be- 
tween A and k. The non-perturbative RG equations in Wegner-Houghton's formulation 

13 



follow if one infinitesimally lowers the running cut-off k — > k — Ak thus retaining the contri- 
bution of the modes which have been previously integrated in the infinitesimal momentum 
shell. In practice the dependent mode improved approximation is obtained by taking the 
derivatives with respect to k of Eq. ( |32| ) rewritten for a general potential and wavefunction 
renormalization. One obtains the following set of non-linear partial differential equations 
for the evolution of Z and U (for simplicity the subscript k is omitted) 

dU lDl Zk 2 + U" - (R/12M2Z + QZ') 

k—— = —ank In- 



dk (Zk 2 + U")\$ =R -_ 



o 



dZ a D k D f9 „ 2Z'K' , Z{K>) 2 , ZZ'k 2 K' Z 2 k 2 (K')\ 

k dk = -—( z - — + ^i- + — 2 K^~) (36) 

where a D = l/2 D n D ^ 2 T(D/2) and 

K = Zk 2 + U" - (R/12)(2Z + <S>Z') (37) 

These equations rule the evolution of all coupling constants generated by the renormalization 
procedure, and have been derived without any assumption on the functional form of the 
potential and the wavefunction renormalization. For instance, if U and Z are analytic 
functionals we write 

OO 1 OO 1 

Z=Y. —h tJ (k)W, U=Y: —lluikWIV. (38) 



k 



gs in 


Eq. (H) are 


dg^ 


d i+j U 


dk 






d i+j Z 


dk 


" d^diK vac 



and one obtains an infinite dynamical system that can be conveniently truncated to some 
order to obtain an approximate solution of Eq. fl56|) . Let us for instance neglect the field 
dependence in Z using the improved version of Eqs. ( p3[ ) in D = 4 (the notation Z = rjf. is 
adopted analogously to Eqs. (|33| ) ) 
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k d£ = __k*_ ^ k k 2 + U" - (R/6)rj k 
dk " 16vr 2 n ( % P + f/")|$= R =o 

d\n Vk _ A; 4 (C/'") 2 (C/" - (i2/6)%) 



(40) 



dink 16tt 2 (% A; 2 + U" - (#/6)?7fc) 4 
It is interesting to employ these equations in order to study the behavior of the couplings 
in the model (|34] ) which corresponds to a specific truncation in Eq. (p8|). 

We define the running parameters m| = #20? = #40, = #21, f]k = h 00 and, as usual, 
the renormalized parameters rriR, Xr, £r, 77^ are obtained in the limit k — > and, in the 
broken phase where (<fi(x)) = v, m 2 R = <?2o( K = 0) + Xrv 2 /2. 



By fixing the vacuum in Eq. (|39|) at the values (<j>(x)) = v and the Ricci scalar at the 



value R = 1Z we obtain the following /3-functions 
arc 1D7T 2 

7fc = ThTfc = -To^ /fc4A ^ 2(m ^ + fe ~ %/6) ^ + W/2) (41) 

where 



] k x = Vk k 2 + m 2 k + X k v 2 /2 + (& - Vk/6)n (42) 

From the above equations it is evident that the presence of a nonvanishing expectation value 
of the scalar field v has the consequence of generating new couplings r k for the operator $ 3 
and a k for $>R. In principle the location of the vacuum should be determined after solving 
the equation of the minimum for the effective potential, in the improved scheme. This 
amounts to solve the system ( fill) coupled to the equation of the minimum for the running 
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potential. Here things are more complicate since in higher-derivative gravitational theories 
the vacuum does not necessarely coincide with the 1Z = configuration • m t ne following 
we just consider v and 71 as adjustable parameters. 

We now discuss the scaling behavior of the couplings in the two extreme scaling limits: 
ultraviolet and infrared. In the UV domain we neglect all scales with respect to k, obtaining 



k 



drrii k 2 X 



k 



dk ?7fcl67r 

r 

k 



d\k 3A| 



dk rjllQ-K 2 

d£ k _ A fc fo - 7fc/6) 

dk t]116ti 2 

dr k 3A^f 

dk r]ll6ii 2 

dak _ vXk f& Vk- 

dk ?7fcl67r 2 



k 
k 

UU k VA k , Ilk, / y| Q \ 

k — = ZJT7^^--) (43) 



while the anomalous dimension vanishes as 

Ik ~ O(p) (44) 

as it should be. This result coincides with the findings in [f7||13l for the /3-functions in the UV 
domain. However once k becomes much smaller than m k , the scaling laws change. Suppose 
first that the following hierarchy holds 

m 2 R ^>7Z (45) 

then, below the mass gap the renormalized flow is arrested and in k = the flow always 
reaches a (completely trivial) fixed point. In particular this shows that the findings in [19| 



that £r = 1 /6 is not necessarely an IR attractor, is still valid when the contribution of rj k is 
considered. At this level of the approximation we note that ^ = 1/6 is a fixed point only if 
?7fe = 1, statement which is not true in the IR domain for the broken phase. In the symmetric 
phase, instead if we fix r/k at the bare level to one it does not change and ^ = 1/6 is a fixed 
point for the system fl4ip . However, as we shall see in the next section, when the contribution 
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of the irrelevant operators is included, this fixed point disappears. These conclusions should 
apply for any D as one can see by direct calculation of the /3-functions derived from Eq. 
(f!0|). We also note that the gaussian fixed point is always present in any dimension. 

When Eq.p5|) is not satisfied the structure of the vacuum can be different. In fact in 
general the curvature competes to determine the sign of the gravitational running mass d\U 
and it can lead the sistem to a new phase. In fact, by looking at Eq. (42) we note that 
around 

Vkk 2 ~ (Vk/Q - £k)n - m% (46) 

the equations in (|41]) become unstable because of the presence of poles in the /3-functions 
signaling that a new phase sets and that new IR relevant interactions are needed to describe 
the renormalized system. 

One might worry that around k 2 ~ R the local momentum expansion is not reliable. 



This is not the case. In fact for the Einstein Universe we have |?7J] the following exact 
finite-difference RG equation for the local potential 

k 

U n -An.k-Ak - U n>h = -~r^n 2 \n(a 2 k 2 + a 2 dlU(<S>, a 2 ) + n 2 - 1) AnAfc (47) 

where n is the quantum number associated with the smearing on the spatial sections with 
topology S 3 , k is a cut-off in the Euclidean "time" direction and a 2 = 6/R is the radius of 
S 3 . We find also in this case the appearance of a singular behavior at some scale n if the 
"restoring force" d^U(§, a 2 ) becomes negative. 

The physical reason of this instability is that the saddle-point expansion in Eq. ([TT|) has 
been performed by considering perturbations around homogeneus configurations which do 
not dominate the path-integral in Eq. flTT| ) on scales smaller than the symmetry-breaking 
scale. In this case a more refined computation is needed, and the naive gradient expansion 
cannot be applied. 
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IV. THE CROSSOVER 



The study of the crossover is instead more complicated and one has to integrate numer- 
ically the flow equations by keeping the contribution of the "irrelevant" operators as well. 
In order to deal with a more tractable problem, we employ a truncation in the expressions 
of the local potential and the wavefunction renormalization function, by keeping up to 
dimension-6 operators. We thus consider the following action defined at the scale k 

S = J d D x^g{e + e,R + |i? 2 + gi? 3 + ^ Vk d^d^ + -^ 2 d^d^ 

+ ^ k Rd^d^ + X -m\& + l -i k m 2 + iA fc $ 4 + ^XkR 2 ® 2 + |i?$ 4 + ^ fc $ 6 } (48) 

where eo = m 2 PL Acs/8n, ei = — mp L /167r, mpi is the Planck mass, Acs is the cosmological 
constant. Note that the remark at the end of sect. II also applies for the more general action 

In order to avoid problems with triviality we keep the UV cut-off A fixed. It is convenient 
to work in running cut-off units, by introducing the dimensionless variables 



x 



^ D - 2 ^ 2 , y = Rk 2 - D (49) 



and calling u, z respectively the dimensionless potential and wave function renormalization, 
we have the expansions 

z = E 7^ z dk)^y j , « = £ -^Uij{k)xy (50) 

i,3 i,3 

where one can easily find the proper correspondence among %(fc), Uij(k) and the various 



couplings appearing in Eq. (E3). If we define 



A t (t, x, y) = z(t, x, y) + d 2 x u{t, x, y) - y{2z + xd x z)/12 (51) 

the renormalization group equations for u and z can be rewritten in terms of the new 
variables in the following way 
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du _{D- 2) 
~dl ~ 2 



(xd x u + 2yd y u) — Du — a D \ia( 



A, 



x=y=0 



dz _(D- 2) 
~dl ~ 2 



(xd x z + 2yd y z) 



a D A~ t \d 2 x z - 2A; 1 d x zd x A t + Afz[{d x A t f + d x zd x A t ] - Af z 2 (d x A t ) 2 ) . (52) 



We have rescaled the momentum variable, by defining t = ln(k/A) so that we generate the 
renormalized flow by lowering the cut-off A — > Ae*, with t < 0. 



$ = at non-zero curvature R = TZ (the corresponding dimensionless value of y is indicated 
with r: r = 1Zk 2 ~ D ), obtaining in D=4 the system of ordinary differential equations dis- 
played in the appendix. For simplicity we shall not discuss here the symmetry broken phase 
with d^O. 

The results of the numerical analysis is shown in figures 1-5, for some values (including 
zero) of the curvature r. We set all bare non-gaussian couplings to zero. The bare value 
of Afc and rjk are kept fixed respectively to u^oit = 0) = 0.1 and zoo(t = 0) = 1, while we 
have considered several bare values of the conformal coupling. We locate the critical line for 
r = near U2o{t = 0) = —0.000325 and this value does not depend on the value of U2\- 

When the theory is away from the critical line, one sees that below a certain value of 
t, the evolution of all couplings stops. In fact below that scale f t , defined in the appendix, 
exponentially vanishes and the tree level scaling of the couplings is recovered. In particular 
for the mass 



The specific constant value of the renormalized mass rriR depends on how the other param- 
eters are fixed. For instance, for the above initial conditions we see that the mass scale 
corresponds to t ~ 6. 

Let us consider the scaling above the mass threshold. The running of some non-gaussian 
couplings is characterized by the presence of a "plateau" as it is shown in fig.l for the running 



We insert Eq. flSTp in Eq. ( |52"D and we evaluate the derivatives in the symmetric vacuum 



; 2 2 




(53) 
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of U03 = e^k 2 and we also observe a similar behavior for U41 = (kk 2 . This is interesting, we 
believe, because it contrasts with the usual notion of "irrelevance" since those interactions 
behave like marginals for almost three orders of magnitude below the cutoff. 

A striking result concerns the running of the u 22 = Xkk 2 interaction term. In fig. 2 we see 
that across and above the mass scale its value is not negligible, since for instance Xkk 2 ~ 0.2 
around t = —4. Note the strong dependence of Xk from the bare value of 

The conformal coupling = u 2 \ deserves a more detailed discussion. Its evolution 
is shown in fig. 3 for the minimally coupled bare theory and in fig. 4 for the conformally 
coupled bare theory. It is worth to remark the small increase of U21 in fig. 3, even for non 
vanishing non-renormalizable couplings (see curve 4), within the three orders of magnitude 
range between t — and the mass threshold, below which the curves become flat. This 
means that starting at t — with a bare £ far from the conformal value 1/6, it is very likely 
that the coupling does not reach that value, being stopped before by the mass threshold. 

Actually, looking at the evolution equation for u 2 \ in the appendix, we see that in the 
infrared region, where the various scales cannot be neglected when compared to k, it is 
extremely difficult to establish the existence of a non trivial IR fixed point for this coupling. 

In fig. 4 the behavior of u 21 in the neighbourhood of 1/6 is explored. Here we see that 
there is no deviation from 1/6 only for zero non-renormalizable couplings at t — 0. However, 
turning on these couplings, we observe small deviations, of the same order of magnitude of 
the changes observed in fig. 3. Again no strong attraction toward u 2i = 1/6 is observed. 

Finally in fig. 5 the influence of the curvature on the phase diagram for the $ field 
is investigated. For a small value of the curvature r = 0.01, we see that for sufficiently 
high values of the conformal coupling (greater than about 2 in our example) the running 
adimensional mass u 2 o(t) grows negative in the IR region while choosing smaller or negative 
u 2 i(t = 0), we see that u 20 (t) is driven away from the critical line toward positive values, 
at earlier "times" t. In fact, for such curvature, the mass in fig.5 practically represents the 
second derivative of the potential with respect to the scalar field, since the additional term 
(£fc — Vk/Q) r is practically t-independent and negligible for large values of t. Thus we are 
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led to a result which is opposite to the naive statement, obtained by looking at the not 
RG-improved propagator, that the phase transition, when r is positive, is approached by 
decreasing the conformal coupling. 

V. SUMMARY AND CONCLUSIONS 

In this paper we have discussed the renormalization of a scalar theory coupled to the 
gravitational field by means of the blocking procedure. We have obtained non-perturbative 
flow equations for the blocked potential and the wavefunction renormalization that rule 
the flow of all coupling constants generated when the cut-off is lowered from the UV region 
towards the infrared domain. We have seen that the gravitational field influences the scaling 
laws of the field in the crossover region between the mass gap and the cut-off. In fact our 
equations reproduce the standard perturbative /^-functions only in the UV region when the 
contribution of the irrelevant couplings is neglected. However we show that any coupling of 
the kind R l & is generated in both potential and wavefunction renormalization even if it is 
not present at the bare level. We have seen that they can be important in determining the 
scaling law across the mass gap. In particular, due to these non renormalizable couplings 
one cannot claim that £ = 1/6 is a fixed point for the conformal coupling constant whereas, 
independently on the bare value of only very small changes of the value of this coupling 
are observed in the flow from the UV to the IR region. 

Another interesting issue addressed with our formalism concerns the structure of the 
phase diagram for the quantum field. In the model analyzed in the previous section the 
onset of criticality is determined by taking into account the coupling with the gravitational 
field. It turns out that the approach to the critical line strongly depends on the strenght of 
the conformal coupling. 

It would be important to understand the consequences of our findings in the contex 
of the inflationary models. For instance it has been recently pointed out that the "fine 
tuning" problem in models with the symmetry-breaking scale near the Planck mass, can be 
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avoided in potentials dominated by non-quadratic "irrelevant" interactions terms like m 
with m > 2. 

In fact in the standard slow-rolling approximation it is assumed that the inflaton field 
evolves as a free field leading to a gaussian spectrum for the large scale angular anisotropy. 
Since this happens above the mass scale of the inflaton field, one may think that the pres- 
ence of non-quadratic interaction terms whose fluctuations are not suppressed above the 
mass scale generates non-gaussian fluctations in the CMBR. At the present state of the 
experimental data there is no secure bound on non-gaussian fluctuations on large angular 
scales, therefore this possibility cannot be ruled out. It would be worth to pursue these 
investigations in a specific model in order to answer these questions in detail. 
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VI. APPENDIX 

In this appendix the complete set of equations ruling the couplings of the action (|48l), 
is displayed in terms of dimensionless quantities introduced in Eqs. ( f49j^0|) for the vacuum 
characterized by the values (4>(x)) = and R = TZ or, in terms of dimensionless quantities, 
x = and y — r. Defining f t as 

ft 1 = u 22 r 2 /8 + u 20 + («2i - z 00 /6)r + z Q0 + z 01 r(l - r/6) (54) 

we get 
~df ~ 2rZ01 ~ 16^ ttZ2 ° 

dz2o 5 1 

— ^ = 2^20 + T^^2o/ t 2 («4o + u 4 ir + z 2 o(l - r/3)) + ^jfH u w + u 41 r + z 2 o{l - r/3)) • 
at lo7T z an z 
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(z 00 + rz 01 )[f t (zoo + rz i)(u 40 + u 41 r + z 20 (l - r/3)) - ru 41 - u 40 - 2 20 (2 - r/3)] 

^01 1 , 2 ( Z 00 U 22 T /-. /Q u 

~dT = 01 16^* Z20 ^ 21 ~ ~6" ~T~ ^ ~ r ' ^ 

= -2u 20 + u 22 — - -J^(u 40 + u 41 r + z 20 (l - r/3)) 
at 4 \bn z 

du 2 i r f t f? 

~W = U22 2 + 16^ (Z2o/3 ~ M4l) + 16^ (M4 ° + W + 

z 20 (l - r/3))(u 21 - Z -f + + z 01 (l - r/3)) 

du 4 Q ft 3f? , , /o\\2 

—jj- = 2ru 41 - Y^ M eo + Y^i( M4 o + M4ir + ^ov 1 ~ r / 3 )) 

<^ M 60 15 2 



16tt 2 



- 2m 60 + jg^/t «6o(«40 + M 4ir + 2 20 (1 - r/3)) - 
30 

f?(u 40 + u 41 r + (1 - r/3)) 3 



— ^ = m 22 /2 + — ^(m 40 + w 4 ir + ^2o(l - r/3))(w 22 /4 - z 10 /3) + 
F ';(u 21 - Z -^ + ^f + z m {l - r/3))[(t*4i - W3) - 



8tt 2V 6 4 

/t(«4o + w 4 ir + 220 (1 - r/3))(«2i - TT + -y- + ^oi(l - r/3))] 

o 4 

= 2m 4 i + ^( M 40 + «4ir + ^20 (1 - r/3))(u 41 - 2 20 /3) + 

/ 2 
f / ^00 . ^22^ 



1C ~ ~g~ + "X" + ^ 01 ^ ~ r / 3 ^[ Meo ~ 6 ^( M4 o + M 4ir + 220 (1 - r/3)) ] 
= -2u 10 - ^(tfci - T + — + *bi(l - r/3)) 

^ M 02 /t /^10 M 22 / 4 2 o M22r , , ,s 2 

~dT = 2w + I6^ ( T - - } + l6^ (M21 - T + — + Z0l{1 - r/3)) 

du os . ^f? ( zoo , u 22 r 3 

~df = 4M3 ° - I6^ (M21 - T + T + Zm{1 ~ r/3)) " 



16tt 2V 3 4 yv z,± 6 

/I o , M ° 3 3 , 1 i ft 

"ST = ~ 2M ° ir + ^ + 16^ ln / t (x = 0,y = r) - (55) 
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FIGURE CAPTIONS 



Figure 1. uo 3 (t) with u 2i (t = 0) = 1, (1), and u 2 i(t = 0) = —1, (2). For both curves the 
non renormalizable couplings are fixed to zero at t — and r = 0.01. 

Figure 2. u 22 (t) with u 21 (t = 0) = -8 and r = (1), u 21 (t = 0) = -8 and r = 0.01 
(2), u 2 i(t = 0) = — 4 and r = 0.1 (3). The non renormalizable couplings are fixed to zero at 
t = for all curves. 

Figure 3. Flow of u 2 \{t) starting at u 2 \(t = 0) = for zero non-renormalizable couplings 
at t — and r = (1), r = 0.01 (2), r = 0.1 (3); with the various non-renormalizable 
couplings set to 0.01 at t — and r = 0.01 (4). 

Figure 4. Flow of u 2 i(t) starting at u 2 i(t = 0) = 1/6 with the various non-renormalizable 
couplings set to 0.01 at t = and and r = (1), r = 0.01 (3), r = 0.1 (4). The flat curve 
(2) corresponds to zero non-renormalizable couplings at t — and it is not sensitive to r. 

Figure 5. u 20 (t) plotted for three different values of the conformal coupling: u 2 i(t = 0) = 
(1), u 21 (t = 0) = 1 (2), u 21 (t = 0) = 2.5 (3); r is kept fixed to r = 0.01 and the the bare 
non-renormalizable couplings are zero. 
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